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ABSTRACT 

We explore the impact of ultraviolet (UV) radiation from massive Population III (Pop 
III) stars of 25, 40, 80, and 120 M Q on the subsequent Pop III star formation. In 
this paper, particular attention is paid to the dependence of radiative feedback on the 
mass of source Pop III star. UV radiation from the source star can work to impede the 
secondary star formation through the photoheating and photodissociation processes. 
Recently, Susa & Umemura (2006) have shown that the ionizing radiation alleviates 
the negative effect by ^-dissociating radiation from 12OM0 PopIII star, since an H2 
shell formed ahead of an ionizing front can effectively shield HVdissociating radia- 
tion. On the other hand, it is expected that the negative feedback by HVdissociating 
radiation can be predominant if a source star is less massive, since a ratio of the 
H2-dissociating photon number to the ionizing photon number becomes higher. In 
order to investigate the radiative feedback effects from such less massive stars, we 
perform three-dimensional radiation hydrodynamic simulations, incorporating the ra- 
diative transfer effect of ionizing and ^-dissociating radiation. As a result, we find 
that if a source star is less massive than as 25M Q , the ionizing radiation cannot sup- 
press the negative feedback of HVdissociating radiation. Therefore, the fate of the 
neighboring clouds around such less massive stars is determined solely by the flux of 
H2-dissociating radiation from source stars. With making analytic estimates of H2 shell 
formation and its shielding effect, we derive the criteria for radiation hydrodynamic 
feedback depending on the source star mass. 
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1 INTRODUCTION 

The reionization and metal enrichment of the universe are 
thought to be gin with the formation of first metal-free 
(Pop III) stars (lGnedinll2000l ; ICiardi et al.ll200ll ; ICenll2003l ; 
ISokasian et alj 2004). Hence, the formation rate of Pop 
III stars is crucial for the subsequent structure formation 
in the universe. The Pop III objects are expected to col- 
lapse at 20 < z < 30, forming a min ihalo with a mass of 
rts 10 6 M B and an extent of sa lOOpc dTegmark et alJll997t 
iFuller fc Couchmanl2000l ; lYoshida et al.ll2003l ). In the course 
of bottom-up structure formation, such Pop III minihaloes 
merge to form first galaxies at z > 10, having the virial tem- 
perature Tvir > 10 4 K and the mass > 10 8 M Q . Even in the 
evolution of first galaxies, Pop III stars can play a significant 
role, since an appreciable number of stars may form from 
meta l -free component in interstellar gas dTornatore et al.l 
120071 ; I Johnson et alJlioog ). 

The formation of very first stars has been investi- 
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gated intensively in the last decade. Many studies have 
come to a similar conclusion that such stars form in a 
top-heavy mass function with th e peak of ~ lOOMp) (e.g., 



Abel. Bryan fc Norman! |2000| ; [5 romm. Coppi fc Larson! 

20021 ; iNakamura fc Umemura! [2001I ; lYoshida et alj 12006). 
Recently. lO'Shea fc Normanl (|2007l ) have shown that the 
variations of cosmological density fluctuations allow the 
mass of Pop III stars to be down to ~ 2OM0 . 

On the other hand, the secondary Pop III star formation 
has been investigated recently. The formation of secondary 
stars is subject to various feedback effects by first stars. One 
of them is the supernova (SN) feedback through mechani- 
cal and chemical effects. The negative feedback by SNe is 
the evaporation of neighboring clouds, since the SN shock 
heats up the gas in clouds. On the hand, SNe can bring 
positive feedback through the compression by shock and the 
cooling by ejected heavy elements. The secondary star for- 
mation can be p r omoted by such positive feedback effects 
jMori et al.ll2002l ; iBromm et alj[2003t iKitavama fc Yoshidal 
120051 ; I Greif et al. 1 120071 ). Another important feedback effect 
is brought by the ultraviolet (UV) radiation from first stars, 
since they are very luminous at ultraviolet band. First stars 
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photoionize and photoheat the surrounding media, and also 
photodissociate H2 molecules, which are the main coolant 
of primordial gas. The radiative feedback from first stars is 
the primary feedback until first stars end the life-time of 
~ 10V with SNe. 

The photodissociation of H2 molecules leads to a neg- 
ative radiative feedba ck effect, which has bee n studied by 
many authors so far. lOmukai fc Nishil (|l999h investigated 
the effect of Eb-dissociating radiation from a single Pop 
III star residing in a virialized halo. They found that if 
the halo is uniform, H2 molecules in the halo are totally 
dissociated, so t h at the gas cannot collapse to form stars. 
iGlover fc Brandl <| 200 if ) considered more realistic clumpy 
halos. They found that if the gas density is sufficiently 
high, photodissociation process proceeds slower than the 
collapse of the cloud. Hence, the cloud can form stars. 
This result is confirmed b y the recent 3D radiation hy- 
drodynamic simulation by ISusal (|2007T ) including the ef- 
fects of hydrodynamics as well as the radiation trans- 
fer of H2-dissociating radiation. The feedback effects by 
diffuse H2-dissociating radiation ca n be important after 
the local feedback in minihaloes jjHaiman, Rees fc Loebl 



Table 1. Properties of Pop III source stars 



19971; iMachacek, Bryan fc Abell l200ll ; lYoshida et all 120031 ; 



O'Shea fc Normanll200^) 



These works basically focused on the photodissociation 
effects. We also have to take into account the effects of ion- 
izing photons. Ionizing radiation heats up the gas through 
the photoionization processes. The temperature of photo- 
heated gas is kept to be around 10 4 K, owing to the bal- 
ance between the radiative cooling and photoheating. If 
the gravitational potential of star forming halos are not 
so deep as to retain the photo heated gas, the heated gas 
evaporates from the halos (e.g., Susa fc Umemural 2004a 3 
Yoshida et all l2007bl ; I Wise fc Abel I 120081 ; IWhalen et all 



2008). However, the case in which ionizing radiation is cou- 



pled with H2-dissociating radiation is complex. When an 
ionization front (I-front) propagates in a collapsing core, 
the enhanced fraction of electrons catalyzes H 2 formation 
JShapiro fc Kandll987l; iKang fc Shapiro I Il992l ; ISusa et all 
ll998l ; IOh fc Haimanll2002r i. In particular, the mild ionization 
ahead of the I-front generates an H2 shell, which potentially 
shield s H2 dissociating photons (|Ricotti. Gnedin. fc Shulll 
120011 ). This mechanism is likely to work positively to 
form Pop III stars. On the other hand, the I-front can 
be accompanied w i th a shock for an optically-thick cloud 
|Susa fc Umemural [2006). The shock affects significantly 
the collapse of cloud. This is a totally radiation hydrody- 
namic (RHD) process. Such radiation hydrodynamic feed- 
back ha s been investigated b y ID spherical RHD simu- 
lations (lAhn fc Shapiro! 120071 ). 2D cylindrical RHD sim- 
ulations (IWhalen et al. [2008), and 3D RHD simulations 
l|Susa fc Umemural 12009 ). The results by 2D and 3D sim- 
ulations are in good agreement with each other. It is found 
that ionizing radiation can bring positive feedback through 
the formation of H2 shell. 

I Susa fc Umemural (|2006t ) investigated RHD feedback 
by a 120Mq source star, and ISusa. Umemura fc Hasegawal 
(2008) derived the feedback criterion. However, if a source 
star is less massive, the relative intensity of H2-dissociating 
radiation to ionizing radiation increases. Then, the feedback 
tends to be more negative. In fact, the mass of first stars 
might be some 1OM owing to the variations of cosmo- 



Mass 


T cS [K] 


JVion [S- 1 ] 


£<Lw[erg/s] 


120M Q 


9.57 x 10 4 


1.069 x 10 50 


5.34 x 10 23 


8OM 


9.33 x 10 4 


5.938 x 10 49 


3.05 x 10 23 


4OM 


7.94 x 10 4 


1.873 x 10 49 


1.17 x 10 23 


25M 


7.08 x 10 4 


5.446 x 10 4s 


3.94 x 10 22 



logical density fluctuations (jO'Shea fc NormanI |2007 | ). the 
enha n ced H2 cooling i n pre-ionized gas JShapiro fc Kand 
1 19871 ; ISusa et al.l Il998l ; lOh fc Haimanl I2002T). or the HP 
cooling in fossil HII regions (e.g . . lUehara fc Inutsuka 
|2000| ; iNakamura fc Umemural |2002| ; Nagakura & Omukai 



2003; I Johnson fc Bromml |2006|; iGreif fc Bromml I200I 
lYoshida. Omukai fc Hernauistll2007al ). Also, the element al 
abundance patterns of hyper-metal-poor stars well match 
the yields b y supernova explosions with a progenitor mas s 
of ~ 25M (jUmeda fc Nomotdl2003l : llwamoto et al.ll2005l ). 
The RHD feedback effects by Pop III stars less massive than 
1OOM have not been investigated so far, and no criterion 
has not been derived. 

In this paper, we perform 3D RHD simulations in or- 
der to investigate the radiative feedback effects from Pop III 
stars with various masses. We derive the criteria for the col- 
lapse of cloud cores irradiated by a neighboring Pop III star 
with 25,40,80, or 12OM . In §2, the simulation code and pro- 
cedure are described. The simulation results are presented 
in §3. Finally, we summarize the conclusions in §4. 



2 SIMULATION CODE AND PROCEDURE 

We perform RHD simulations with a 3D Radiation-SPH 
code developed by ourselves. In the code, we treat self- 
consistently the gravitational force, hydrodynamics, the ra- 
diative transfer of UV photons, non-equilibrium chemistry 
for e", H+, H H", H 2 , and Rj . W e use the chemical net- 
work solver in lKitavama e t al. (200 l|) a s well as the radiative 
transfer solver described in ISusa! (|2006). For the shielding by 
H2 molecules against H2-dissociating radiation at Lyman- 
Werner (LW) band (11.26-13.6 eV), we employ the self - 
shielding function introduced by iDraine fc Bertoldil (|1996| ). 
The opacity against LW band flux (Flw) is calculated by 



-FlW = Fuw,ofs (A r H 2 ,14) 



(1) 



where _Flw,o is the incident flux, A r H 2 ,i4 is the H2 column 
density in units of 10 14 cm -2 , and 



-3/4 



x s; 1 

X > 1 



(2) 



In this paper, we simulate the evolution of a purely 
baryonic primord ia l clou d, according with the model by 
ISusa fc Umemural (2006). The cloud is initially uniform 
with the density of hh = 14cm -3 , and has the mass of 
M = 8.3 x 1O 4 M . The initial chemical compositions are 
assumed to be the cosmological compositions provided by 
iGalli fc Paiial j|l998l ). Before the UV irradiation, the cloud 
contracts self-gravitationally to form a collapsing core. As 
for the core temperature T c , we employ two models. One 
is a high temperature model, and the other is a low tem- 
perature model. By changing the initial temperature of the 
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Figure 1. Evolution of the core without radiative feedback. The 
time variations of core temperature are shown for the initial tem- 
perature 350K (red curve) and 100K (blue curve), respectively. 



clouds Xini, we realize such core temperatures. As shown in 
Fig. [T] if we set T[ n i — 100K, the core temperature becomes 
T c ~ 300-400K at core density n c > 10 2 cm" 3 . On the other 
hand, if the initial temperature is set to be 2I n j = 350K, the 
cloud core cools below T c ~ 200K, since H2 molecules are 
rapidly formed owing to the high initial temperature (see 
FigHJ. Another difference between two models is the ratio 
of gravitational energy W to internal energy U , because it 
is dependent on the initial temperature. The [W|/{7 ratio is 
~ 4 for the low initial temperature (high T c ) model, while 
the |W|/C/ ratio is « 2 for the high initial temperature (low 
T c ) model. 

We ignite a source star when the density of cloud core 
exceeds a certain value n on . The source star is placed D pc 
away from the center of cloud core. We change the mass of 
source star in the range of 25A^0 ^ M* 120Mq . The prop- 
erties of source stars as the effective temperature of star T e g , 
the number of ionizing photons emitted per second A i Qn , an d 
the luminosity at LW band are taken from lSchaererl |2002), 
which are summarized in Table 1. Note that we do not con- 
sider the lifetimes of source stars in this paper, since we focus 
on elucidating the RHD feedback before SN explosions. 

Numerical runs are characterized by the parameters D, 
n on , and M*. The simulations are performed until t en d = 
2tg, where iff is the free-fall time determined by n on . If the 
density of cloud core exceeds 5 x 10 5 cm -3 before t cn d, we 
stop the calculation, since the cloud is expected to keep col- 
lapsing. In order to clarify the effects of the ionizing radia- 
tion, we also perform the simulations artificially disregard- 
ing ionizing radiation but still including LW radiation, and 
the results are compared with those of normal simulations. 
The number of SPH particles handled in our simulations is 
262,144 for all runs. 

The present simulations are mainly carried out with a 
novel hybrid computer system in University of Tsukuba, 
called FIRST simulator, which has been designed to sim- 
ulate multi-component self-gravitat ing radiation hydrody - 
namic systems with high accuracy (|Umemura et al.ll2007t ). 
The FIRST simulator is composed of 256 nodes with dual 
Xeon processors, and each node possesses a Blade-GRAPE 
board, on which GRAPE-6 chips, that is, the accelerator 
of gravity calculations, are implemented. The peak perfor- 
mance of FIRST simulator is 36.1 Tfiops. 
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Figure 2. Time variations of density profiles along the sym- 
metry axis for a simulation with T c ~ 400K, M» = 8OM0, 
n on = 10 3 cm — 3 , and D = 40pc. The red, green, and blue lines 
represent the profiles at lMyr, 2Myr, and 3 Myr after the ignition, 
respectively. The results without ionizing radiation are shown by 
dashed lines. It is shown that the core cannot collapse without 
ionizing radiation. 



1 Myr after ignition 
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Figure 3. Various physical quantities along the symmetry axis 
at lMyr after the ignition for the simulation shown in Fig. f3] The 
results with and without ionizing radiation are shown by solid 
and dashed lines, respectively. The magenta, blue, red, green, and 
black lines show the gas temperature(K), number density (cm -3 ), 
electron fraction, H2 column density from the source star in units 
of 10 14 cm -2 , and the H2 fraction, respectively. 

3 TYPICAL RESULTS 

In this section, we show the typical evolution of clouds. 
For a high core temperature model (T c ~ 300 — 400K), 
the time evolution of density profiles along the symmetry 
axis is shown in Fig. [2j where the set-up parameters are 
M* = 8OM , n on = 10 3 cm~ 3 , and D = 40pc. In this figure, 
the results with ionizing radiation are compared to those 
without ionizing radiation. In the simulation with ionizing 
radiation, the density of cloud core keeps increasing, and the 
density exceeds the limit (5 x 10 5 cm~ 3 ) due to the run-away 
collapse at 3.4Myr after the ignition of the source star. On 
the other hand, in the simulation without ionizing radiation, 
the gravitational contraction of cloud core is stopped by the 
thermal pressure, and eventually a hydrostatic core forms. 

Various physical quantities along the symmetry axis at 
lMyr are shown in Fig. [3] If ionizing radiation is included, a 
dense shell forms ahead of the ionization front (I-front) . The 
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1 Myr after ignition 
2Myr after ignition 
3Myr after ignition 




10 10 
Distance [pc] 

Figure 4. Same as Fig. [2] but for M, = 25M . D = 14pc is set 
up so that the LW band flux toward the cloud core should be the 
same as Fig. [2] 



1 Myr after ignition 
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Figure 5. Same as Fig. [3] except that the parameters are M* 
25M Q and D = 14pc. 



H2 molecule fraction is raised up to a level yH 2 ~ 10 -5 in the 
shell, so that H2 column density exceeds 10 14 cm -2 . Owing 
to the self-shielding of LW band radiation by the shell, the 
H2 fraction in the cloud core is increased, compared to the 
case without ionizing radiation. Eventually, the enhanced 
H2 cooling allows the core to undergo the run-away collapse. 
On the other hand, unless ionizing radiation is included, LW 
band radiation from the source star reduces the H2 fraction, 
so that the cloud core is settled in a hydrostatic configura- 
tio n. This mechanism is ba si cally the same as that found 
by ISusa fc Umemural l|2006| ); ISusa. Umemura fc Hasegawal 
(|2008h in the case of M* = 12OM . 

In Fig. [4] and Fig. [5] the results in the case of M 4 = 
25Mq are shown, where n on = 10 3 cm -3 . Here, the source 
distance is set to be D = 14pc so that the LW band flux 
toward the cloud core should be the same as that in the 
case of M» = 80Mq, whereas the flux of ionizing radiation 
is about 0.75 times weaker than that in the M* = 8OAf0 
case. As shown in Fig. [4] the cloud fails to collapse and 
a hydrostatic core forms, notwithstanding the presence of 
ionizing radiation. Similar to the case of M, — 8OM0, the 
H2 fraction ahead of the I-front is enhanced associated with 
a dense shell. However, the H2 column density of the shell is 
not high enough to shield the H2 dissociating photons. Thus, 
the H2 fraction at the cloud core stays as low as yu 2 ~ 10 -6 , 



10- 



10 J 



E 



10 s 



10' 



10 



— 1 Myr after ignition 




— 2Myr after ignition 




. — 3Myr after ignition 











40 30 20 10 10 20 30 40 
Distance [pc] 

Figure 6. Same as Fig. [2] but for T c ~ 200K. 
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Figure 7. Same as Fig. [3] but for T c ~ 200K. 

which is the almost same level in the case without ionizing 
radiation. Thus, in this lower stellar mass case, the ionizing 
radiation cannot suppress the negative feedback. 

In Fig. [7J and Fig. [7] the results for the low core temper- 
ature model (T c ~ 200K) are shown. In these simulations, 
the parameters are set to be M* = 80M©, D — 40pc, and 
n on = 10 3 cm -3 . As shown in Fig. [S] the cloud fails to col- 
lapse, despite the presence of ionizing radiation. However, 
the reason for the failure is different from the M* = 25Mq 
case. It can be seen in Fig. [7] that the shielding effect raises 
the H2 column density, compared to the case with no ion- 
izing radiation. In this low core temperature model, hydro- 
gen molecules are strongly destroyed by the LW radiation, 
and the H2 fraction decreases to yH 2 ~ W~ 8 . Since the 
core radius is smaller for the low core temperature, the self- 
shielding for LW radiation by the core is weaker. 



4 CRITERIA FOR RADIATIVE FEEDBACK 
4.1 Numerical Criteria 

In Fig. [8] the numerical results are summarized for a high 
core temperature model (T c ~ 300 — 400K). In this fig- 
ure, crosses denote the failed collapse, triangles represent 
the successful collapse with the aid of ionizing radiation, 
and circles represent the collapse regardless of ionizing radi- 
ation. As shown in this figure, in the simulation runs with 
M* > 4OM0, the H2 shell driven by ionizing radiation can 
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Figure 8. Numerical results for T c ~ 400K are summarized in terms of D and n on . From the top-left panel to the bottom-right panel, 
each panel represents the result with M* = 120Mq, 80Mq, 40Mq, and 25Mq. Crosses denote the runs in which the clouds fail to 
collapse, triangles represent the runs in which the clouds can collapse if the ionizing radiation is included, and circles represent the runs 
in which the clouds can collapse even if the ionizing radiation is not included. Thick solid lines indicate the analytic criteria D cr .d> which 
is derived from balancing photodissociation timescale to free-fall timescale (see the text for the detail). On the other hand, dotted lines 
indicate the analytic criteria D CI s j, including the shielding effect by the shell. In each panel, 0.4D cr ^ and 0.4D cr s j, are also shown by a 
thin solid line and a thin dotted line, respectively. 



allow the clouds to collapse if the conditions for D and n on 
are satisfied. However, in the case of M, = 25Mq, ionizing 
radiation does not help the clouds to collapse, but the fate 
of clouds is determined solely by Hh-dissociating radiation. 
Hence, we conclude that the critical stellar mass below which 
ionizing radiation cannot extinguish the negative feedback 
by photodissociation is M* ~ 25Mq. 

In Fig. [9] the numerical results are summarized for a low 
core temperature model (T c ~ 200K). The tendency is qual- 
itatively the same as the results for T c ~ 300 — 400K. But, 
the regions of the collapse with the aid of ionizing radiation 
(triangles) are obviously narrower, and in wider regions the 
clouds fail to collapse. These results basically originates in 
the fact that the radius of cloud core is smaller, compared 
to a high core temperature model with T c ~ 300 — 400K. 
For the smaller core radius, H2-dissociating radiation is li- 
able to permeate and suppress the core collapse. However, 
the critical stellar mass, M* ~ 25M@, below which ionizing 
radiation cannot extinguish the negative feedback by pho- 
todissociation, is almost the same as that in the case with 
T c ~ 400K. This fact means that the critical stellar mass 
does not depend sensitively on the cloud core temperature. 
As a result, we conclude that the critical stellar mass below 
which ionizing radiation is not important is M* ~ 25Mq. 



4.2 Analytic Criteria 

Here , we m ake analytic estimation of the feedback criteria. 
ISusal (|2007l ) explored the photodissociation feedback of a 



Pop III star with 120Mq on a neighboring prestellar core by 
RHP simul ations which dose not include ionizing radiation. 
ISusal (|2007h has found that a condition for the collapse of 
a neighboring core is approximately determined by tdis = 
tg, where fdi s is the photodissociation timescale in the core 
and is is the free-fall timescale. Using the condition, the 
critical distance -D cr ,d, below which a neighboring core fails 
to collapse, is given by 

r> — U7 I Llw i 2 ( nc \ ~ 15 
cr ' d ~ P %5xlO»»ergB-V VlO^m^ 

K (5^c)"*' (3) 

where Llw, n c , and T c are the LW luminosity of source 
star, the number density of core, and the temperature of 
core, respectively. This equation involves the self-shielding 
effect by the core. The dependence on the core temperature 
basically originates in the core radius (oc T^ 2 ) and the H2 
formation rate in the core (oc T c ). Hence, the self-shielding 
effect is wea ker for the lower core temperature. As argued in 
ISusal (|2007l ). the boundary between the collapses regardless 
of ionizing radiation (circles) and with the aid of ionizing 
radiation (triangles) in Fig. [8] can be roughly explained by 
-D cr ,d in the case with M* = 120-Mq. In addition, as shown 
in Figs. [8] and [9] D CI> ^ gives a good estimate for less mas- 
sive source star cases. However, as shown in Figs. [8] the 
boundary for the high core temperature model is slightly 
lower than this analytic estimate. This disagreement can be 
understand by the dynamical effect of the collapsing clouds 
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Figure 9. Numerical results for T c ~ 200K are summarized in terms of D and n on . From the top-left panel to the bottom-right panel, 
each panel represents the result with M* = 120Mq, 80Mq, 40Mq, and 25Mq. Symbols have the same meanings as those in Fig. [8] 



l|Susall2007l ). The actual dynamical contraction is faster for 
the high core temperature model, since the ratio of grav- 
itational energy to internal energy is higher (|W|/!7 « 4) 
as described in Sj2] Then, the H2 fraction in the core recov- 
ers rapidly, during the adiabatic compression phase. Hence, 
the core can keep collapsing, even if the photodissociation 
timescale is shorter than the free-fall timescale (i d is < iff) 
when the cloud irradiated by UV. As a result, the criterion 
is reduced to /dynAr.d by a dynamical factor /dyn- Com- 
pared to the numerical results, we find /d yn ~ 0.4 for the 
high core temperature model. On the other hand, for the 
low core temperature model, the dynamical effect is not so 
strong because of |W|/(7 ~ 2, and therefore /d yn ~ 1. 

Furthermore, if the ionizing radiation is included, we 
should incorporate the shielding effect by an H2 shell. Here, 
we derive a new criterion including the this effect. Since a 
cloud collapses in a self-similar fashion before UV irradia- 
tion, the density profile of outer envelope in the cloud is 
expressed as 



(4) 



where r c is the core radius which roughly corresponds to the 
Jeans scale; 



1 / 7rfcsT c 



Grrili 



(5) 



where fee denotes the Boltzmann constant and m p denotes 
the proton mass. Assuming that the thickness of the H2 shell 
is determined by the amount of ionized gas in the envelope 
and the H2 fraction in the shell is constant, the H2 column 



density of the shell A r H 2 ,sh is given by 



■V, 



H 2 ,sh — 



yn 2 , s hn(r)dr = yH 2 , s hn c r' ( 



DD sh 



(6) 



where -D s h, and yn 2 ,sh are the distance between the cloud 
core and the H2 shell, and the H2 fraction in the shell, re- 
spectively. Here, D s h is set to be the position where the 
number of recombination per unit time in the ionized region 
around a source star balances with the number rate of inci- 
dent ionizing photons, since the H2 shell appears ahead of 
ionization front. Hence, D s h satisfies 



iVion7r£) s 2 h 

4tt(D - Ah) 



27to:b 



n(r) r dr 



2 4 D — L> s h 
— 2naBn c r c 



(7) 

where as is the recombination coefficient to all excited levels 
of hydrogen. Using equation © and J7} , we obtain 



N11 



12 _2 / JV io 



\87TQE 



Because of the intense LW radiation, the H2 at the shell is 
in chemical equilibrium. Therefore, yu 2 ,sh is given by 



2/h 2 



n(D B h)y c ,shk H - 



(9) 



where y e ,sh is the electron fraction at the H2 shell and fe H - is 
the reaction rate of H~ process. In this case, we should con- 
sider the self-shielding effect by the shell itself. As a result, 
these rates are 



1.0 x 10" 18 T sh cm" 3 s _1 , 



(10) 



Radiative Regulation of Population III Star Formation 7 



fc dis = 1.13 x 10 8 F LW , sh / s 



jV H2 , s h/2 
10 14 cm- 2 



(11) 



where T s h and .FLw.sh are the temperature at the shell, and 
the LW flux from the star in the absence of shielding effect, 
-FLw.sh = £lw/47t(D — Ash) 2 - In addition, / s is the self- 
shielding function given by ©. Combining equations @- 
(|11|) with assumption of y e ,sh = 0.05 and T sn = 2000K as 
shown in the present numerical results, we have 



2/H 2 ,sh 



1.0 x 10~ 



300K 



10 50 s -l 

i ( D 



5 x 10 23 erg s~ 



V40pc 



(12) 



H 2 ,sh 



= 5.8 x 10 



v 5 x 10 23 erg S" 1 7 ' ^ 

Notice that Nn 2 ,sh is determined solely by the ratio of Ni on 
to Z/lw, and strongly depends on the ratio. In the above nu- 
merical results, it is shown that the critical stellar mass does 
not depend sensitively on the cloud core temperature. This 
fact is consistent with equation (|13p . in which the H2 col- 
umn density of shell is independent of the core temperature 
T . 

Multiplying L LW in equation © by / SjSh = 
f s ( loi^cm-^ ) ' we °bl am l ne critical distance as 



A: 



= 147pc 



Lhwfs, 



5 x 10 23 erg s _1 



/ n c \ - Te 
V10 3 cm- 3 / 



\300K J 



(14) 



in which both shielding effects by the core and the H2 
shell are taken into account. In particular, when A r H 2 , s h > 



10 i4 cnT 



D B 



the critical distance can be expressed as 
Luw \ 2 / A'ion 



78.8pc 



Vl0 3 cm- 3 / 



5 x 10 23 erg s _1 

16 / 2c 



10 5C V 



V300K 



(15) 



In Fig. [5] and Fig. [5] triangles are the collapse with the aid 
of an H2 shell. Therefore, the boundary between the trian- 
gles and the crosses should be compared with 73 C r,sh- Ac- 
cording to equation (|13p . the shielding effect by the shell 
becomes weaker according as Ni on /Lmv decreases. This in- 
dicates that -Dcr.sh approaches D C r,d as the mass of source 
star becomes lower. As shown in Fig.JS] -D cr ,sh is the almost 
same as D C r,d in the case with M* = 25Mq. This result 
originate in the strong dependence of the shell H2 column 
density on Mon/.I'LW (see equation 1 13p . As shown in Figs. [5] 
and[9l -D cr ,sh gives a qualitatively good estimate for the col- 
lapse with the aid of an H2 shell. However, the boundary in 
Figs. [H] is slightly lower than this analytic estimate. For the 
same reason as in the case of Ar.d, the dynamical effect is 
more prominent for the high core temperature model. In this 
case, /d yn Ar,sh with / dyll ~ 0.4 provides a more appropriate 
criterion. On the other hand, for the low core temperature 
model, /d yn ~ 1 gives a plausible criterion. 



5 CONCLUSIONS AND DISCUSSION 

We have carried out RHD simulations to investigate the im- 
pact of UV radiation from a Pop III star on nearby collaps- 
ing cores. In particular, our attention has been paid to the 
dependence of UV feedback on the mass of Pop III star. 
The radiation hydrodynamic evolution of cloud core is de- 
termined by not only liVdissociating radiation but also ion- 
izing radiation. As a result, we have found the critical stel- 
lar mass M* « 25M , above which an H2 shell ahead of 
ionizing front can help clouds to collapse. Owing to the fact 
that Eb-dissociating radiation becomes predominant for less 
massive source stars, the critical distance for the collapse 
of a neighboring core does not so strongly depend on the 
mass of source star. Also, we have derived analytically the 
feedback criterion, fdynD cr ,sh, where -D cr ,sh is given by (|14|l 
and / dyll is a dynamical factor dependent on the the ratio 
of gravitational energy W to internal energy U of collaps- 
ing cloud. We have found /d yn ~ 0.4 for |W|/{7 « 4, and 
/dyn ~ 1 for \W\/U w 2. Since /d yn is dependent on |W|/f , 
a dark matter (DM) halo can influence the feedback crite- 
rion to a certain degree. In order to assess the effects of 
DM, we have calculated several m odels with a static NFW- 
type dark matter halo potential (|Navarro. Frenk fc White] 
1 1997ft with M vir = 4.15 x 1O 5 M and r vir = 160pc. In these 
runs, the ratios of DM mass (Mum) to baryonic mass (Mb) 
at the central regions of r < lOpc are Mum/M^ ~ 0.3 for 
n on = 10 3 cnT 3 , and M OM /M h ~ 1 for n on = 10 2 cm~ 3 . As 
a result, we have found that the feedback criterion in the 
form of /dynAr.sh turns out to be still valid, and /d yn be- 
comes smaller by a factor of 1.2 for n on = 10 3 cm -3 and by 
a factor of 2 for n on = 10 2 cm -3 . Therefore, our main re- 
sults are not changed so much by including DM. Note that 
the DM density evolution is not treated consistently with 
the gas dynamics in these simulations. If the DM dynamics 
is solved with the evolution of gas clouds, the evolutionary 
path of core temperature might be changed. Hence, for a 
more quantitative argument, the self-consistent treatment 
of dark matter would be requisite. 

In this paper, we have not considered the lifetime of 
source stars. The lifetime of Pop III star is 2.5 x 10 6 yr for 
12OM , 3.0 x 10 6 yr for 80Mp, 3.9 x 10 6 yr for 4OM , and 
6.5xlO 6 yrfor25M (|Schaerer!l2002l ). If the lifetime of source 
star is shorter than the free-fall time determined by n on , 
the feedback may be significantly changed before the cloud 
collapse. The density in which the free-fall time equals the 
stellar lifetime is n on = 419cm~ 3 for 12OAf , n on = 293cm -3 
for 80M Q) n on = 178cm" 3 for 4OM , and n on = 64cm" 3 
for 25M . Below these densities, arguments including the 
effects from the stellar lifetime are requisite. 

The fate of Pop HI stars depends on the mass 



epe 

jHeeer fc Wooslevl |2002| ; iHeger et al.ll2003ft . Pop III stars 
with 12OM or 8OM may result in direct collapse to black 
holes (BHs), while those with 40 M or 25M may undergo 
Type II supernova explosions. In the case of direct BH for- 
mation, UV source disappears abruptly, and then already- 
formed H2 molecules can promote the collapse of cloud 



fe.g.. [Nagakura fc Omukail |2005 | ; Johnson fc Bromm 



20061; jGreif fc Bromml 120061 ; lYoshida. Omukai fc Hernquist 
2007aft . In the case of Type II SN explosions, shock- 
driven hydrodynamic feedbacks could be significant 



8 K. Hasegawa, M. Umemura and H. Susa 



dMori et al.ll2002|; iBromm et all 120031 ; iKitavama fc Yoshidal 
120051 ; iGreif et alj|2007h ~ 
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